On the size of the smallest scales in cosmic string networks 
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We present a method for the calculation of the gravitational back reaction cutoff on the smallest 
scales of cosmic string networks taking into account that not all modes on strings interact with all 
other modes. This results in a small scale structure cutoff that is sensitive to the initial spectrum 
of perturbations present on strings. From a simple model, we compute the cutoffs in radiation- and 
matter-dominated universes. 
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I. INTRODUCTION 
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Topological defects are a generic prediction of grand 
unified theories. They can be formed in phase transitions 
when the topology of the vacuum manifold of the low 
energy theory is non-trivial M . For a review see O . 

The result of a phase transition that produces strings is 
a network of long strings that stretch across the horizon 
and a collection of closed loops. If the phase transition 
occurs at energy r; s , then the mass per unit length in the 
strings is /i ~ rQ. Immediately after formation, strings 
in the early universe undergo an epoch of heavy damping 
until a time £*, which depends on the string formation 
scale rj s as well as the specifics of the particle physics 
model that produces them. As a result of this damping 
all sub-horizon features, such as loops and perturbations 
on long strings, are wiped out and the long strings and 
loops are Brownian with a persistence length d ~ f » . 

After that time, numerical simulations show that in 
an expanding universe there is an attractor solution in 
which the network evolves into a "scaling regime" (see 
|^| and references therein), where the energy density of 
the string network is a small constant fraction of the ra- 
diation or matter density and the statistical properties of 
the system, such as the correlation lengths of long strings 
and average sizes of loops, scale with the cosmic time t. 
This solution is possible because of intercommutations 
that produce cosmic string loops which in turn decay by 
radiating gravitationally. Simulations also have found 
that most loops and the perturbations on long strings 
have the smallest possible size, the simulation resolution, 
which does not scale. The prevailing opinion on this is- 
sue is that the size of small-scale structure in fact also 
scales with the cosmic time t and its value is given by 
the gravitational back-reaction scale TGfit, where T is a 
number of order 100 and G is Newton's constant. This 
possibility was first pointed out in H . 

Cosmic strings are good candidates for a variety of 
interesting cosmological phenomena such as gamma ray 
bursts , gravitational waves §,§ and ultr a high energy 
cosmic rays Some of the predictions of these mod- 



els, however, depend sensitively on the so far unresolved 
question of the size of the small-scale structure. 



II. BACK-REACTION MODEL 

Hindmarsh |J (see also |l(J) showed that the power 
per unit length radiated into gravitational waves from 
two colliding arbitrary small perturbations of long repeat 
length L on an infinite string is 
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where K n — 2im/L and n m — 2irm/ L are the wavenum- 
bers of the right- and left-moving Fourier modes that 
make up the arbitrary perturbations and e n and e m their 
(small) amplitude to wavelength ratios. 

We can use this expression to construct a simple back- 
reaction model. If we split the sum of the wavenumbers 
in Eq. © 
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we can identify individual modes with the power they 
radiate 
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with a similar expression for the m modes travelling in 
the opposite direction. Each of the modes contributes 
to the effective mass per unit length of the string by an 
amount 
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and the power radiated into gravitational waves by each 
of the modes decreases this contribution [9| 
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By putting Eq. (||) together with Eqs. (Q) and (||) we 
arrive at an expression for the time evolution of the am- 
plitude to wavelength ratio of each of the modes 
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with an analogous expression for the left-moving m 
modes. Each mode therefore loses amplitude at a rate 
proportional to its own frequency, to its amplitude to 
wavelength ratio and to a sum which represents the in- 
teraction of the mode with every other mode moving in 
the opposite direction. 

There remains, however, some puzzling behavior. In 
the case of just two colliding modes we can write Eq. ([!]) 
as 
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where the a and b subscripts differentiate between the 
right- and left-moving modes. If we conformally stretch 
(say) the a mode such that e a remains constant while 
K a —> 0, Eq. (0) approaches the constant 
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This result is in contradiction to the fact that the power 
vanishes in the case of perturbations travelling in only 
one direction on a string [ fi"l| : By conformally stretching 
one of the modes we are making the string on which 
the other mode is travelling straighter and we expect the 
power radiated to approach zero in that limit. 

This problem was investigated in jl2|]. Here we will 
summarise the results. It turns out that in order to be in 
the Hindmarsh regime (the regime where Eq. (Q) is valid) 
it is not merely sufficient for the amplitudes of each of 
the colliding modes to be smaller than their correspond- 
ing wavelengths: It is also necessary for both amplitudes 
to be small compared to the geometric mean of the wave- 
lengths, namely it is necessary to have 
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where A a and Ab are the amplitudes and A a and A& are 
the wavelengths of the two colliding modes. If instead we 
are in the regime where one of these two Lorentz invariant 
quantities, say, 
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the power radiated is exponentially suppressed. This is 
precisely the regime where the power given by Eq. (Q) 



tends to the constant Eq. (g). It was also found that 
when the quantity e 2 Kb/SK a increases past 1 the power 
drops discontinuously. 

We can approximate this behaviour by introducing a 
cutoff in Eq. (Q) such that the sum excludes interactions 
between modes that satisfy 
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as follows 
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with a similar equation for the m-modes traveling in the 
opposite direction. This approximation is justified be- 
cause of the discontinuous drop in the power that takes 
place when e 2 Kb/K a = 8 and the exponential suppres- 
sion that sets in when Eq. ( p"l| ) is satisfied. In fact we 
cannot easily associate the power radiated with a given 
mode outside of the Hindmarsh regime and it is fortu- 
nate that in this case the power radiated is exponentially 
supressed and therefore negligible. It should be noted 
that the calculation in |lj] was performed for just two 
modes travelling in opposite directions. Here we have 
assumed (reasonably, we believe) that when a mode in- 
teracts with a collection of modes all of which satisfy Eq. 
( fil"| ) the power is also negligible. 

Eq. ( |13| ) can be used to find the evolution of a set of 
amplitude to wavelength ratios from some initial condi- 
tions. There are two distinct cases of interest. On a long 
string, one can consider Eq. ([l3]) to give the damping 
of the average perturbation with a certain wavelength. 
By symmetry the left-moving and right-moving pertur- 
bations will have the same average amplitudes. However, 
in the case of a loop, a statistical fluctuation may pro- 
duce an excess of power in a particular direction over a 
broad range of wavelengths. In that case, since a pertur- 
bation passes the same oppositely-directed perturbations 
over and over again, it is essential to treat the right- and 
left-movers separately. Here we will consider the case of 
long strings only. 



III. THE SIZE OF THE SMALLEST SCALES 

When the time-scale over which modes on a long string 
are significantly stretched by the expansion of the uni- 
verse, which is about t, the cosmic time, is large com- 
pared to the gravitational damping time, Eq. ( |T3| ) can 
be used to estimate the size of the smallest scales on long 
strings: The nth mode will be significantly damped in 
amplitude provided r g < t with r g given by 
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We expect that a combination of stretching by the ex- 
pansion of the universe and self-intersections on strings 
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will result in an average spectrum of pertubations where 
the small scales are suppressed relative to the large scales. 
This seems reasonable even when kinks (produced at in- 
tercommutation events) are present since the amplitude 
to wavelength ratio of Fourier components of kinks as a 
function of mode number decays like 1/n. The presence 
of a flat or decreasing spectrum means the sum in Eq. 
( |l3| ) is going to be bounded from below and above: A 
given mode does not see every other mode on the string 
but rather a range of modes that lie around it. If we 
assume a power law spectrum 

e„~n- /3 , (15) 

one can use Eq. ([l2]) to show that the lower bound on m 
for a given mode n is 
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which corresponds to the largest wavelength mode that 
the nth mode interacts with while still in the Hindmarsh 
regime. The smallest wavelength mode that it interacts 
with can also be found from Eq. (n% and is given by 



m Max ~ 8n 1+2 ^. 
We can therefore write Eq. (|13f) as 
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with mMin and mMax given by Eqs. fllq ) and (|17|) above. 

Previously it was assumed that all modes interact with 
all other modes and therefore for spectra where the small 
scales are sufficiently suppressed, namely those with (3 > 
1/2, 
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This means that we would expect the survival of modes 
whose wavelengths 
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independently of the specific spectrum of perturbations 
present on the string. 

In fact we need to evaluate the sum between the two 
limits. We can approximate the sum in Eq. ( [l8|) by an 
integral 
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and take the dominant contribution which when (3 > 1/2 
is given by 
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The dominant contribution therefore comes from the 
largest wavelength mode that the nth mode can inter- 
act with. 

Using Eq. (|22|) we can write Eq. (O) as 
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Using Eqs. ( |l4| ) and (16) we can see that 
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Since we want the gravitational lifetime of the mode to 
be larger than the cosmic time and using n ~ t/X, 
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and therefore only modes with wavelengths 

A>(rc A1 ) (1+2 ' 3)/2 i, 

have a significant amplitude at time t. 
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IV. SIMPLE EXAMPLES 

When the amplitude is small compared to the wave- 
length, one can find the effect of cosmological expansion 
on the spectrum of perturbations on a string Jig] ] . If the 
wavelength of a mode on a long string is larger than the 
horizon, then it is stretched conformally; its amplitude 
and wavelength grow with the scale factor. If, on the 
other hand, the wavelength of the mode is smaller then 
the horizon, only its wavelength is stretched. The net re- 
sult of this cosmological processing is that the amplitude 
to wavelength ratio of modes larger than the horizon as 
a function of the wavenumber remains constant and for 
modes inside the horizon it becomes a power law. 

It is not hard to see why this is so. If a mode enters the 
horizon at a time to, then its wavelength and amplitude 
are both ~ to. While inside the horizon the amplitude 
remains fixed and the wavelength is stretched so that at 
sometime t later 
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where a — 1/2 or 2/3 depending on whether we are in 
the radiation or matter era. We can therefore write the 
amplitude to wavelength ratio of that mode as 



eo(t) 



• n 



t]T a t a 



(29) 



3 



If we imagine that our mode of wavelength Ao (t) is some 
fraction of the largest mode, Ao(i) — X/n with A ~ t, it 
is easy to see that 

t Q = tn- 1 ^ 1 -^ (30) 
which when substituted into Eq. (^9|) yields the spectrum 

_1 Radiation Era 
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yields 



Feeding the spectrum Eq. (31) into Eq. 

A > (TG )( 1+Q )/ 2 ( 1_Q ) t — < O^m) * Radiation Era 
~ \ (TGfi) 5/2 t Matter Era 

(32) 

as the small scale structure cutoffs in the radiation and 
matter eras. However, this analysis neglects intercom- 
mutations, which produce kinks, and the possibility that 
small perturbations propagating on curved horizon-sized 
strings may not be efficiently stretched. 



V. CONCLUSIONS 

We have presented a method for the calculation of 
gravitational back-reaction on cosmic strings taking into 
account that not all Fourier modes that make up the per- 
turbations on strings interact with all other modes with 
the same efficiency Q . 

In particular, modes of a given wavelength only in- 
teract significantly with a narrow range of other modes 
whose wavelengths are comparable to it. This range de- 
pends on the spectrum of perturbations on the string at 
late times (when gravitational effects become important) 
which in turn has the effect of making the small scale 
structure cutoff sensitive to this spectrum. 

Using this method and assuming a power law spectrum 
we have arrived at an expression for the small scale struc- 
ture cutoff. We have further applied it to the spectrum 
resulting from the stretching of small amplitude waves in 
a radiation- or matter-dominated FRW universe. 

It is unclear to what extent the results obtained for 
these two simple cases apply to a realistic network of 
cosmic strings because we have ignored the effect of in- 
tercommutations, which leads to kinks, as well as the pos- 
sibility that small perturbations propagating on curved 
horizon-sized strings may not be significantly stretched 
by the expansion of the universe. Both of these effects 
have the potential to change the spectrum of perturba- 
tions and therefore the value of the small scale structure 
cutoff. 

However, it is clear that because of the restricted range 
of interaction affecting every mode, the small scale struc- 
ture cutoff will be smaller than the one given by the usual 
calculation where an unrestricted range of interaction is 
assumed. 
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